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QQ • Abstract. We give a complete analysis of low- velocity dynamics close to grazing 

for a generic one degree of freedom impact oscillator. This includes nondegenerate 
(quadratic) grazing and minimally degenerate (cubic) grazing, corresponding 
Q . respectively to nondegenerate and degenerate chatter. We also describe the dynamics 

r \ I associated with generic one-parameter bifurcation at a more degenerate (quartic) graze, 

showing in particular how this gives rise to the often-observed highly convoluted 
structure in the stable manifolds of chattering orbits. The approach adopted is 
geometric, using methods from singularity theory. 
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O . 1. Introduction 

The theory of forced nonhnear oscillations has well known and widespread applications 

/\ ' in science and engineering: see Hayashi [2lj, Fidlin |20], for example, as well as landmark 
^^ ' 

cd I texts Andronov et al. 0, Schmidt et al. [37], Guckenheimer et al. [23]. However, 

there are many engineering contexts in which oscillations are a highly undesirable 
feature of performance, especially when these lead to impacts between moving parts or 
against rigid supports, generating noise and wear and ultimately leading to mechanical 
breakdown: for some examples see Alzate et al. p]. Mason et al. [31j, Stewart |40j . 
Pavlovskaia et al. [35] as well as those described in books such as Babitsky et al. [5], 
Brogliato [Hj and di Bernardo et al. [7]. 

A simple model for such dynamical impact phenomena (often called vibro-impact 
systems) is the one degree of freedom impact oscillator. This consists of a forced 
nonlinear oscillator 

x + f{x,x)=g{t), g{t + T)=g{t) (1) 
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for X G R and fixed T > 0, toEretfier with a constraint x > c where c is a constant often 
called the clearance. Here / and ,, are smootlfl toction^T and a= u=ual the dot denotes 

d/dt. 

The model must include a prescription for how the dynamics are to proceed when 
X = c. This will typically be given by a restitution rule that whenever (x, x) = (c, v) 
then the velocity v is instantaneously replaced by v (a function of v and perhaps other 
variables) satisfying vv < 0. There are many variations on this model, involving more 
than one degree of freedom and/or more than one constraint or a non-instantaneous 
restitution rule. 

The fact that a simple impact system can give rise to extremely complicated 
dynamics has long been noted, and there are many analytical and numerical studies 
showing complicated periodicities or convoluted orbit structures in phase space: see 
for example Thompson et al. [44J, Shaw et al. [38], Hindmarsh et al. [25], Peterka et 
al. [36], Kleczka et al. [30j, Virgin et al. L46j, Wagg et al. [1Z],[18], Valente et al. [15]. 
Nevertheless, even a 1-degree of freedom linear oscillator with a simple restitution rule 
is still poorly understood from the overall dynamical point of view. 

A more systematic attack on the global problem was undertaken by Whiston in a 
pioneering series of papers [19] , [50] , [HI] ; this enterprise inspired much subsequent work 
including Budd et al. [15], Lamba [28], Nordmark [33] as well as the present paper. 

A key feature of impact oscillator dynamics that gives it a different flavour from 
(smooth) nonlinear dynamics is the phenomenon of grazing where a periodic orbit 
reaches the obstacle with zero velocity and nonzero acceleration: small perturbations of 
such orbits may or may not exhibit impacts. The resulting bifurcation phenomena have 
been studied by many authors including Budd et al. [H], Chin et al. [TTj, Dankowitz et 
al. [18], Foale et al. [21], Ivanov [2Sj, Nordmark [33], [3l], Szalai et al. [12], Molenaar et 
al. [32], Zhao et al. [53j . 

Another important aspect of impact oscillator dynamics that attracts particular 
interest is the low- velocity behaviour close to impact, characterised by the phenomenon 
of chatter. On leaving x = c with small positive velocity but negative acceleration the 
orbit (i.e. the solution to ([T])) may return to hit x = c again after a short time; after 
restitution at the impact the process will repeat, and under suitable conditions the orbit 
may come to rest at x = c after a finite time but an infinite number of such impacts. 
It will remain stuck at the obstacle x = c until such time as the acceleration (from the 
forcing) becomes positive and it once again lifts off into the region x > c. See for example 
Demeio et al. [19] for a simple mechanical realisation of this phenomenon, or Stone et 
al. [H] , Valente et al. [15] for more complicated dynamical settings. The first systematic 
analytic and geometrical analysis of the dynamics of chatter was carried out by Budd 
et al. [H], using a combination of insightful numerics supported theoretically by first 



X For simplicity we take smooth to mean C°° , although many of the results will clearly hold with less 
differentiability. 
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and second order approximations. In this paper we extend those results by providing a 
complete and rigorous analysis of chatter for generic impact oscilaltor systems ([1]) and 
for generic 1-parameter families of such systems. 

A standard technique for studying a forced oscillator such as ([1]) without the impacts 
is to focus on the the discrete dynamical system F that samples the solutions at times 
{nT : n G Z}: then periodic solutions of ([T]) correspond to periodic orbits of F, and 
so on. When impacts are introduced, however, this is no longer appropriate since the 
influence of the obstacle on a given trajectory cannot be expected to be T-periodic. A 
common approach to dealing with this problem, used by Shaw et al. [38], Thompson 
et al. [B], Whiston [50] and others, is first to regard (^ as an autonomous system 
in R^ by including the equation t = 1 and then to consider the plane x = c and 
the discrete dynamical system defined on it by the first return of trajectories after 
impacting the obstacle. This fits naturally with the geometry of the physical description 
of ([1]), but the resulting discrete system now has loci of discontinuity that arise from 
grazing trajectories. It is not easy in this setting to see how the discontinuities and 
the associated dynamical phenomena respond to change in initial conditions and/or 
governing parameters such as the clearance. 

In this paper we use an alternative geometric approach, first set out in [16] . 
which allows greater insight into the geometry of discontinuities and their influence on 
dynamics by replacing the picture of curved orbits intersecting a flat obstacle with an 
equivalent picture of straightened-out orbits intersecting a curved obstacle: the impact 
surface. See Section [2] below for the precise description. From this point of view, the 
problem of analysing the geometry of orbits at and near grazing becomes more readily 
amenable to the tools of local differential topology, dynamical systems and singularity 
theory. 

Taking this approach, we look again at the phenomena of grazing and chatter in 
the light of the geometry of the impact surface. We flrst use this to provide rigorous 
underpinning for the key results in [13] that are based on approximations: by placing 
the results in a wider context we are able to give a full analysis of the local dynamical 
phenomena. We then extend this to cover the case of the discontinuous bifurcations 
that arise in the neighbourhood of a doubly-degenerate graze as the clearance or some 
other parameter is varied and the degeneracy is unfolded. In particular we exhibit 
the key mechanism for the creation of multiple-loop structures in stable manifolds of 
chatter points, which is a flrst step in understanding the multiple-loop patterns observed 
numerically for attractors and basins of attraction of periodic orbits in impact oscillators: 
see [501 Figure 17], [HI Figures 4 and 19], as well as [H Figure 3] and [7] (in particular 
the front cover). Taken together, the results of the present paper provide a complete 
description of generic low-velocity impact dynamics and 1-parameter bifurcations. 
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2. Geometrical model for dynamics: the impact surface 

We begin by recalling the construction described in Chillingworth pL6j for modelling 
the overall dynamics of a one degree of freedom impact oscillator. The model is based 
on a 2-dimensional surface in three dimensions that characterizes the impacts, together 
with certain natural maps associated with it that represent the dynamics. A similar 
geometrical construction was used by Sotomayor et al. [39], albeit with a different 
dynamical interpretation, to classify local vector fields on a manifold with boundary. 

2.1. The impact surface and associated maps 

Let x{b, V, r; t) denote the unique solution to ([1]) with initial data 

(x, x) = (6, v) when t = t. 

For fixed choice of c G R we then write 

Xc{r, V] t) := x{c, v,T]t + t) 
so that 

Xc{t,v;0) = c, Xc{t,v;0) =v. (2) 

Definition 2.1 The impact surface Vc is where Xc — c vanishes, that is 
Vc := {(r, V] t) e R^ : Xc{t, v; t) = c}. 

Note that in [T6| the convention was to write (f , r; t) rather than (r, v;t) : we trust no 
confusion will arise from this notation reversal. 

By definition Vc contains the plane t = 0, which we shall denote by 11 and which 
plays a key role in our formulation of the dynamics. It is proved in fTEj Lemma 1] that Vc 
is indeed a surface (smooth 2- manifold), with the exception of the points on the r-axis 
Ho and possibly certain other points where v = 0, characterised precisely in terms of the 
functions / and (? in ([1]) and the clearance c . It is showr|5| in [T6| Proposition 1] (see 
also [14j where it is inherent in the analysis, and see the calculations below) that in a 
neighbourhood of IIq the impact surface Vc consists of two sheets (smooth 2-manifolds) 
intersecting transversely along IIq . One of these is of course the plane 11 while the other 
sheet V^ is represented in a neighbourhood of IIq as the graph of a smooth function 
V = Vc{T,t) given implicitly by 

yciT,Vc;t) = (3) 

when Xc is written in the form 

Xc = c + tydr, V] t) (4) 

for a smooth function yc on a neighbourhood of IIq in R^. 

S The statement r > 2 there is an error: it originally referred to degree of differentiability, and has no 
connection with the restitution coefficient here. 
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VciT,t)\t=Q = -^acir) (5) 



By implicit differentiation of (^ we find 
d_ 
di 

where adr) := Xc(r, 0; 0) is the initial acceleration (again see [E] as well as [lH 
Proposition 1]). 

There are three maps which we use to build the dynamical model for the impact 
oscillator. We now describe each of these in turn. 

2.1.1. The projection pc Consider the projection 

p : R^ -> n : (r, f ; t) t-^ (r, v) 
from R^ to the plane 11 of initial data, and let p^ denote its restriction to V/ : 

p,:=p\V^:V^^n. 

Let He be the set of singular points of pc : K' — )■ 11 (including those, if any, where V^ 
itself may be singular). From the Implicit Function Theorem it is immediate that 

H, = {{T,v;t) e V; : x,{T,v;t) = 0} 

which in geometric terms is the horizon of V^ as viewed along the t-direction. 
Dynamically, the points of He represent graze points, that is points where the solution 
to (II]) meets the obstacle x = c with zero velocity. See Figure [H where for diagrammatic 
simplicity the curve H^ is not distinguished from the horizon of V^ as seen from the 
viewpoint of the reader: in reality these are of course different. The subset Pc := p{H^ 
of n is the apparent outline (classically the apparent contour) of the surface V^ in the 
t-direction: see Bruce [TOl .Bruce et al. [T2l. 



The generic local geometry of smooth maps from a smooth surface M into R^ is 
well understood since the work of Whitney [52] . The set of singular points is a smooth 
1-manifold in M consisting of arcs of fold points (where the singular set maps locally to 
a smooth 1-manifold) and isolated cusp points (where the singular set maps locally to 
a |-power cusp). The precise characterisations are given in terms of partial drivatives, 
and local coordinates can be chosen so that at fold points or cusp points the map has 
the form (x,y) h-> {x,y'^) or {x,y) h-)- {x,xy + y^) respectively. These generic results 
apply also to the conceivably more restricted setting of projection of a surface in R'^ to 
a plane, and thus to the geometry of apparent outlines. In Figure [1] the points Wi,W3 
are fold points, while W2 is a cusp point. 

As the surface or the view direction varies with one or more parameters the outline 
will typically undergo certain local transitions. The generic 1-parameter transitions in 
apparent outlines are also well understood: see [3|,[9],[T3] and Figure [2 

In general, without any restriction on t, there may be regions of 11 filled densely by 
Pc . However, if Vc is periodic in t (as can certainly happen for linear systems) or we 
restrict to a compact subset of Vc then this complication will not arise. 
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Figure 1. The impact surface Vc — HUV^ with horizon He that is the set of singular 
points of the projection p\V^ — > 11. Points Wi,W2,'W3 belong to the apparent outline 
Pc — p{Hc) (not shown). 

It is pertinent to ask whether for arbitrary choice of clearance c the apparent 
outhne of the impact surface V^ , which is determined by solutions of the differential 
equation ([T]), is indeed generic as just described, and whether the transitions in the 
apparent outline as c varies are also generic in this context. These questions may be 
posed both for a general system ([T]) and for a specific system such as the forced linear 
oscillator 

X + X = COS cut. (6) 

A detailed discussion of the generic geometry of V^ and its outline for the general 
system ([1]) and for the linear system iQ can be found in [16j where these questions are 
answered. 



2.1.2. The re-set map (pc There is a natural smooth re-set map 0c : K — ?• 11 obtained 
simply by re-setting the clock to t = at the moment of impact: 

(pc-Vc^U: (r, v; t) ^ {t + t, Xc{t, v; t)). 

Let Zc denote the zero set of Vc, that is the intersection of V^ with the plane v = 0. The 
following facts about the singularity structure of 0c are proved in p^, Propositions 6,7]. 
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Figure 2. The three generic transitions in apparent outlines: (a) lips, (b) beaks, 
(c) swallowtail. 



Proposition 2.2 The singular set of 0c is Z^ \ Hq- Points (r, 0; t) G Zc\Ilo where 
adr) 7^ are fold points of (pc, while points where adr) = but a'^r) ^ are cusp 
points. 

2.1.3. The restitution map Re The restitution map describes how the velocity 
immediately after impact depends on the velocity immediately before impact. Thus 
we write 

n+ = {{t,v) eU:v>0} 
U- = {{t,v) eU:v<0} 

and we suppose that Re : H" — > 11+ is the restriction to 11^ of a map 

Rc:U^U:{t,v)^ (r, pc{v)) 

where pc : R- —> R- is a differentiable function with Pc(0) = and Pc(0) < 0. In many 
cases an appropriate choice is Pc{v) = —rv with < r < 1 independent of c, although 
all that we shall require is that Pc(0) = —r. 



2.2. The first-hit map F^ and the dynamical system Gc 

Since we are assuming the dynamics take place in the region x > c it is only the upper 
half-plane n^j' := 11+ U Ho that is of physical interest. Nevertheless, in much of what 
follows it will be useful first to consider 11 as a whole and then restrict to 11^ later. 
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The projection pc : Vc ^ H has a partial right inverse constructed as follows. 
Choose initial data (t, t>) G 11+ (so x = c and f > at t = 0), then proceed parallel 
to the t-axis in the direction of increasing t until the point Fc{t,v) of first intersection 
{first hit) with Vc is reached: thus 

F,iT,v) = iT,v;t^)eV^ 

where (r, v;t) ^ Vc ioi < t < ti. It is possible that no such ti exists (no future impacts 
occur), but for notational simplicity we disregard such points and write 

F, : n+ -> V^ 

in any case, with the understanding that Fc may not be defined on the whole of 11"'". 

The construction of Fc extends naturally to points (r, 0) G IIo for which Xc{t, 0; t) > 
for all sufficiently small t > 0; if this does not hold (for example, if Xc{t, 0; 0) < 0) 
then the constrained orbit remains stationary at x = c and we define Fc{t,0) = (t, 0). 
The points of V^ in the image of Fc : 11^ — )• V^ will be called visible points, as these are 
the points of V^ that can be 'seen'from II^j" along the positive t direction. 

Finally, although it is not neded for dynamical purposes, for geometric convenience 
we define Fc : 11" — ?■ V^ analogously, proceeding here in the direction of negative t. 

The first-hit map F^ : 11 — )■ V^ defined in this way typically has discontinuities. 
Since Pc is a local diffeomorphism at points z E V^\Hc it follows that the discontinuities 
of Fc must lie in the apparent outhne Pc = p{Hc). Not all points w E Pc are discontinuity 
points of Fc , of course, but only those with Fc{w) G He : thus in Figured] the point Wi 
is a point of discontinuity of Fc while the point w^ is not. 

2. 3. The correspondence principle 

The task of describing the dynamics of the impact oscillator close to grazing is greatly 
helped by exploiting a correspondence, induced naturally by 0c , between geometric 
features of Pc close to IIo and those away from IIq . 

li z E He then Pc{z) G Pc by definition, and since Xc = at 2; we have (pdz) G IIq. 
Likewise ii z E Zc then Pc{z) G IIq by definition, and it is also the case that ipdz) G Pc ■ 
This is because z represents current data (on x = c) for a solution with initial velocity 
f = so, by uniqueness of solutions of ([1]), after re-set the point (t)c{,z) represents initial 
data for a solution x{t) with x{t) = c and with x{t) = at some other time t . Thus 

(t)c{Hc) = IIo = Pc{Zc) and Pc{Hc) = Pc = MZc)- 

In fact by extending this argument we obtain a more general result which we leave as 
an exercise: 

Proposition 2.3 

p;\PcUUo) = (j)-\PcUUo). 

n 
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Furthermore, since re-setting the time clearly does not affect geometric features of the 
dynamics, any features of the outline Pc in a neighbourhood oi w = Pc{z) where z E He 
must have counterparts in a neighbourhood of 0c(-2) € Ho. This correspondence principle 
between P^ and Ho enables us to focus the analysis on points of Ho when studying local 
dynamics close to grazing, and (as we see below) provides a valuable framework for 
understanding the structure of the dynamics itself. A fuller description of the geometry 
associated to the correspondence principle can be found in [T6] . 

2.4- Construction of the dynamics 

We are now ready to construct the (discontinuous) discrete dynamical system on the 
half-plane 11,} that we will use to model the dynamics of the impact oscillator. 

Begin by applying the first-hit map F^ : TIq — )■ V^. This corresponds to following 
a solution of ([1]) as it leaves the obstacle at x = c until it returns to x = c for the first 
time. Next, apply the re-set map 0c : K' ~^ n; this is simply a re-calibration of the 
time variable. Since the initial data lie in 11^ the velocity at first hit must be negative 
or zero, so we have 

0c o i^c : n+ ^ n^ := n- U Ho. 

Finally apply the restitution map Re : Hq — ?■ Hq : 

The iteration of this process reconstructs the impact dynamics of the oscillator as 
a (discontinuous) discrete dynamical system on Hq generated by the map 

G, ■= i?c o 0c o Fc : n+ -^ n+. 

As indicated above, we may regard F^. and Re as maps defined on the whole of IT, and 
hence the same applies to Gc whenever this is convenient. 

The rest of this paper is devoted to describing fully the generic behaviour of 
iterations of Gc close to the r-axis Ho , both for fixed clearance c and as a 1-parameter 
family of dynamical systems when c or another parameter is varied. We take the partial 
results of [H] further by giving a rigorous derivation of the dynamics for nondegenerate 
chatter at regular points of Ho (that is where adr) ^ 0) and degenerate chatter at 
tangency points of Ho (where ac{T) = but the derivative a'^lr) ^ 0). We also consider 
a yet more degenerate situation that generically arises only for certain values of the 
clearance c at so-called swan points of Ho (where adr) = a^^r) = but a^^r) ^ 0), 
and we show how the dynamics unfold and undergo newly-created discontinuities in 
a 1-parameter bifurcation as c passes through such values. The creation of new cusp 
points for pc gives rise to convoluted stable manifolds of nearby chatter points, with 
immediate implications for complicated global dynamical behaviour. Understanding 
the global structure of stable manifolds is one of the main challenges in formulating a 
dynamical theory for non-smooth systems such as impact oscillators. 
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3. Nondegenerate chatter 

Let the clearance c be fixed. To simphfy notation in this section we temporarily write 
V, F, G etc. for V^, F^, Gc etc., although we retain Vc for clarity and understand that p 
denotes p\V^. Using coordinates (r, v) on 11 and (r, t) on V , we identify a point (r*, 0; 0) 
of the r-axis Ho = 11 fl V with both w* = (r*, 0) G 11 and 2;* = (r*, 0) G V . 

Let w^ = (r*, 0) G Ho be a regular point of Ho , so that a{r^) 7^ 0. Since the 'free' 
dynamics takes place in the region x > c , the behaviour of any solution that leaves 
X = c with zero velocity but positive acceleration will not be susceptible to our local 
(low-velocity) analysis close to Ho : therefore we assume a* := a^r^:) < . On leaving 
X = c with small positive velocity the solution to ([T]) will slow down, then possibly 
reverse and hit x = c again after a short time. This is indeed the case as we see 
from (jS]): since z^ ^ H the map p : V ^ His a local diffeomorphism at z^, and therefore 
F is also a local diffeomorphism at w* . From Proposition 12.21 we have that (p : V ^ U 
is a local diffeomorphism at w^, and soG' = i?o(;/)oFis (the restriction to n^j" of) a local 
diffeomorphism of U at w^ . 

3.1. Local linearisation along Ho 

As a first approximation to the dynamics we linearise at the equilibrium point u7=k G Ho. 
Since by definition 

G op = R o (p 

we have on differentiation 

DG{w)Dp{z) = DR{u)D(t){z) 

for z &V' where w = p{z) and u = (f){z), both in 11. Writing r = r,, + a we see from ([5]) 
that the derivative Dp{z^) is 

Dp{z^) : (a, t) H- (a, -ia*t). 

Also, using implicit differentiation of (jH]) to see that ^{z*) = |a* we find (see [161 
Lemma 2]) 

D(f){z^) : (a, t) i-> (a + t, |a*t). 

Hence the derivative of the local diffeomorphism G = i? o o F at w* is given by the 
2x2 matrix 

DGiw.)=(l "M (7) 

where r = — p'(0) is the restitution coefficient (compare [T^). 

Eigenvectors of DG{w^) are (1, 0) with eigenvalue 1 (since the r-axis Ho is fixed) 
and (2,a=i,(l — r)) with eigenvalue r. Therefore, to first order, we have the picture as 
illustrated in Figure [3](a) (see also [13]): orbits of the local linearisation of the discrete 
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Figure 3. Local dynamics of G in the (r, w)-plane 11 close to Ho: nondegenerate 
chatter for a{T) < 0. (a) Local linearisation, (b) local foliation by stable manifolds. 



dynamics G approach Ho along directions with negative slope, and this slope tends to 
zero as a^, — > or as r — )■ 1. 

This does not complete the story, however, as it remains to confirm that this linear 
description also accurately describes the nonlinear behaviour of G on 11 close to Uq. 

A fundamental component in the theory of smooth dynamical systems is the 
analysis of dynamics in the neighbourhood of a normally hyperbolic invariant manifold, 
that is an invariant manifold for which the rate of normal contraction or expansion 
overpowers any contraction or expansion along the manifold itself: main references 
are |26] and |13]. See also [27] for an encyclopaedic treatment. 

In our case, let J be any open interval of the r-axis Ho on which a{T) is negative 
and bounded away from zero. Then G is the identity on J and at each point w & J 
the derivative DG{w) has a contracting eigenspace (since < r < 1) transverse to Ho. 
Hence J is normally hyperbolic for G, and from the standard theory [26] we immediately 
deduce the following description of the local dynamics of G on some neighbourhood U of 
J in n. Recall, however, that for appplications to the impact oscillator we are concerned 
only with lid', ^^^^ i^ where v > 0. 

Theorem 3.1 For each t ^ J there is a smooth G -invariant curve W^ in U G H 
through Wr = (t, 0) G IIq and tangent to the direction (2, (1 — r)a{T)). The curve W^ is 
a local stable manifold for Wr-' 

W^ = {u e U : G"-{u) -^ Wr as n ^ +oo}. 

The curves {ly^j^gj provide a continuous foliation of U, in that there is a 
homeomorphism $ : J x R — )■ f/ taking {r} x R to W^ for each r E J. See Figure\^(h) . 

Thus the linear picture is indeed an accurate representation of the nonlinear dynamics 
close to Ho , although this result yields only that the stable manifold W^ varies 
continuously with r. Can we obtain more, and see that in fact W^ varies smoothly 
with r so that the union of the W^ form a smooth foliation of a neighbourhood of IIq ? 
From a result of Takens [l3] the answer turns out to be positive. 
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Theorem 3.2 // the functions f and g in IjJ^ are C°° then (for suitable U) the 
homeomorphism $ in Theorem \3.1\ may he taken to he a C°° diffeomorphism. 

Proof. This follows from the more general (although local) theorem in [43j concerning 
partially hyperbolic fixed points. The non-resonance conditions which are imposed in 
the more general setting of |43j are in our case vacuously satisfied since < r < 1 and 
there is no expansion or contraction along the manifold J. Extension from the local 
result to the whole of J follows by usual compactness arguments applied to the closure 
of J. D 

Therefore the dynamics of the local linearisation as described in [T3] do indeed give 
an accurate picture of the true (nonlinear) dynamics up to smooth change of coordinates. 



As aij) — 7- the stable manifold W^ approaches IIo with slope tending to zero. 
Our next step is to investigate the local dynamics in this limiting case. 

4. Degenerate chatter 

In this section we examine the nature of chatter close to a point w* = (r^,, 0) where now 
O'd'T*) = 0. For generic choice of the functions / and (/ in ([1]) and for an open dense set 
of choices of the clearance c the T-periodic function ac{T) will have (modulo T) a finite 
set of zeros each with a'^{T) 7^ 0. In particular this is the case for the linear system (|6]), 
since here ac{T) = —c + cosur where T = 27t/uj. However, for certain values of c the 
function ac{T) can typically be expected to undergo Morse transitions (fold bifurcations) 
where adr) = and a'^^r) = while a"{T) 7^ 0. We shall describe the local dynamics 
of the impact oscillator in both these situations. 

We start by calculating the first few terms of the Taylor series of the maps 

Pc: (r,t) ^^ (r,t;e(r,t)) (8) 

and 

0c : (r,t) h^ (r + t,Xc(r,Wc(^,i)) (9) 

at points (r, 0) G Hq . Since from (jl]) we have 

V = Xc{t,v;0) = yc{T,v;0) (10) 

it follows that -^ = 1 and all other derivatives of yc with respect to v and/or r vanish 
on Uq. Using this, by implicit differentiation of ([3]) and writing r = r^, + a we find 

Vein + a,t) = - |at (11) 

-iaVi + i(-ia' + ia6)t2 (12) 

- \a"ah + i(-|a" + ^a'b)at^ (13) 

- |(i«" + §(-!«' + bbW + 0(4) (14) 
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where a = ac(T*), b = -^Xc{t^, 0; 0) and the primes denote derivatives with respect to r 
evaluated at r^, (which in general depend also on c). Hence, if a 7^ then 

Vc{n + a,t) = -lat + 0{2) (15) 

while if a = but a' ^ then 

Vcir, + (r,t) = -|a't(3a + t) + 0(3). (16) 

Finally, if a = a' = but a" ^ then 

Vein + (T,t) = -l4:a"t{6a^ + Aat + t^) + 0(4). (17) 

4-1. Dynamics near a tangency point 

A point (r, 0) G Ho where ac{T) = and a'^{r) 7^ is called a tangency point. In the 
context of Section [3] this can be regarded as a point of Ho where the normal hyperbolicity 
fails in the least degenerate way. 

Suppose that w^ = (r*, 0) G Ho is a tangency point. From ( !T6|) we see that Vc has a 
saddle critical point at {a, t) = (0, 0). The horizon He is tangent to the line 3cr + 2t = 
which projects from V^ by Pc to the parabola v = |ac(T*)cr^ in 11. We consider two cases 
separately, according to the sign of ^^(t*). Again for convenience we drop the suffix c 
from the notation. 

(a) The case a'^r^,) > 

The points of H close to w^ have v > 0. The part of H relevant for us is the arc 
where t > 0, that is a < 0, and so the relevant part of P = p{H) here is a smooth 
arc P* quadratically tangent to Ho at w^ from the direction r < t^ and v > 0. This 
is the discontinuity set for the first-hit map F and hence for G locally: we call it the 
discontinuity arc. See Figure |5l 

Theorem 4.1 Let w^, = (t*,0) be a tangency point with a'ij^) > 0. The map G 
has a unique invariant arc 7 C Hq with an end-point at w^,, this arc being smooth 
and quadratically tangent to Uq at w^, from the side t < t^ and lying between P* and 
Uq. The fixed point w^ is exponentially attracting for GI7, while 7 is weakly normally 
hyperbolic (repelling) for G. Orbits close to but not on 7 either (if below ^) are attracted 
exponentially to Ho in nondegenerate chatter or (if above '-/) lose contact with the obstacle 
on crossing the discontinuity curve P* after finitely many impacts. See Figure 

By weakly normally hyperbolic we mean that at points w G 7 the derivative DG{w) has 
an eigenvalue A 7^ 1 (in this case A > 1) with eigendirection transverse to 7. However, 
A — 7- 1 as If — )■ tu^, along 7. The terms above and below correspond to values of v. 

A proof of the existence of such an invariant curve 7 was outlined by Budd and Dux 
in [H] by restricting to second-order terms. The analysis in [H] is also for a particular 
linear system, although the authors observe that their results apply more generally. 
We provide here a rigorous proof that applies in full generality to nonlinear systems of 
type (HD. 
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Figure 4. The actions of p, and R close to (r*, 0) with o! (t^^ > 0: here r — t* + cr. 
The shaded area in the top left diagram represents the visible part of V', and the other 
shaded areas are its images under the relevant maps as indicated. 



4-1.1. Proof of Theorem 4-1 From (|8]) and flT6|) we see that the projection p^ : V ^ U 
is given in local coordinates by 

p: {a,t)^ {a,-k{3at + t^) + 0{3)) (18) 

where k = ^a'^r^) > 0, and from ([9]) and (TT6l) the reset map : l^' — )■ IT is given by 



: (a, t)^ {a + t, -k{3at + 2t^) + 0(3)). 



(19) 



In Figure m we indicate some local geometry of p and of 0. As already noted, the horizon 
H is given by 3a + 2t = 0, while to first order the zero set Z is given by 3cr + 1 = 0. 
Also, the line a + 1 = is taken by (p to the f-axis cr = 0. 

Now take new coordinates (a, rj) in R^, where 1] = ^a + t. Then 

p : (a, T]) ^ ((7, -kr]^ + ^ka^ + 0(3)). 

Composing with the diffeomorphism 



^ : R2 ^ R2 : (m, w) ^ (m, w - |A;m^) 
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Figure 5. Degenerate chatter at w^. = (t»,0) with a'{Tf,) > 0. Invariant curves below 
the invariant curve 7 represent nondegenerate chatter; invariant curves above 7 carry 
orbits across the discontinuity arc P,. 



we obtain 

ip o p : (a,ri) \-^ (a, —kri^ + 0(3)) 

which exhibits p exphcitly as a fold map. By standard methods of singularity theory (see 
e.g. Arnol'd et al. [1], Golubitsky et al. [22j), after composing with C°° diffeomorphisms 
of R^ in range and domain that are the identity to second order, we may in fact assume 
that 

ijj op: (a,r]) h^ {a,-kri'^) (20) 

without any higher order terms. 

Let g : y — )■ n denote the composition Rocf). In the new (a, rj) coordinates we find 

V^ o g : (a, r/) ^ (r/ - \a, -rk{2rf - 3ar]) - ^k{r] - ^af + 0(3)) (21) 

where we recall that r = — p'(0). 

To say that 7 is an invariant curve for G is equivalent to saying that the curve 
(3 = F{'y) C V has the property that p{f3) = q{f3), which in turn is equivalent to saying 
that P is an invariant curve for the map h := p~^oq whose domain and range are subsets 
of v. Explicitly we have 

h = p-^oq:{a,r])^iri- \<j, v/Qr(a,r7) + 0(3) ) (22) 

where 

Q,(a, r/) := r{2r,^ - 3av) + f (r/ - ^cxf. (23) 
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Our strategy for finding /3 will be to blow up the coordinates at {a, rj) = (0, 0) and then 
exhibit /3 in the new coordinates as the stable manifold of a hyperbolic fixed point of 
G. 

For the blow-up construction we write 

{a,r]) = (s,sm) =: xis,m) 

so that m is the slope of the ray from the origin on which {a, rj) lies, and s measures 
displacement along the ray. Then 

h{s, sm) = ( (m — ^)s, A/s^Q^^(1^7n)~+7J(s^ ) 
and so using (s, m) as coordinates this becomes 

h' : {s,m) H-)- {s',m') 
where 

s' ={m- |)s (24) 

m' = {m- i)- VQr(l, m) + 0{s). (25) 

The domain Dr of h is the subset of R^ corresponding to the values of m for which 
m 7^ I and the expression under the square root is nonnegative. The following is 
straightforward to check. 

Proposition 4.2 The quadratic function m i— )■ Qr{l,m) has zeros at m^,m^ where 

< m_ < i < m+ < I 

and Qr{l,m) > for m G Ej. := (— oo,m_] U [m+, oo). Also m_, m+ — )■ | as r — > 0. D 

Since distinct zeros of a quadratic function are necessarily simple zeros, the Implicit 
Function Theorem yields the following result. 

Corollary 4.3 For sufficiently small e > the domain Dr of h' consists of the region 
Dr{e) = D;{e) U D+(e) where 

D;{e) := {is,m) eR^ : \s\ < e,m < 5_(s)}, (26) 

D+(e) := {(s,m) G R^ : |s| < e,m > 5+(s)} (27) 

and 6± are smooth functions with S±{0) = m± respectively. 

We now seek a fixed point for /i' on x E'^. This occurs where 

m= {m- \)'^^JQr{l,m), 
that is where 

Qr{l,m) = m^{m-\f. (28) 

Solutions of f l28|) are given by m = | and by the solutions of 

2rm={m-\f{m + l). (29) 



Dynamics of an impact oscillator near a degenerate graze 17 

A simple sketch of the graph of the cubic function on the right hand side of (129 p shows 
that this equation has three solutions mi,m2,m3 with 

mi < — I and < m2 < | < "^3- 

We now investigate these four solutions of (!28l) and verify that only one of them is 
relevant here. 

Solution m = I : in the original (a, t) coordinates this corresponds to t = 0, that is the 
r-axis Ho. Clearly this is invariant under G (on one side of (r*, 0)) but it represents 
equilibrium states and is not the arc 7 representing chatter that we seek. In fact 
771 = I implies r = 1 which is the ideal case of perfect restitution in the limit as 
V ^0. 

Solution TTii < — I : the image of this line is given opposite orientations by p and q, so 
that under the composition q op~^ the orientation of the image is reversed. Such a 
curve cannot be a candidiate for 7. Moreover, since the gradient of the line is less 
than —3 it maps under p (or q) to a curve tangent to the r-axis but with f < , 
not relevant for us. 

Solution < 777-2 < I : the image of this line is also given opposite orientations by p 
and q (in particular, points with f > are taken to points with w < 0) so this image 
cannot be a candidate for 7. 

Solution 7773 > I : here when t > both p and q take this line into the region f > 
with the same orientation, so if 7 = p{P) exists as an invariant curve then the 
tangent direction to /3 at the origin must correspond to this ray. 
Note that 7713 increases from | to | as r increases from to 1. 

We next study local linearisation of h' at (0, m^). Using (|22|) and (|23|1 we find that the 
Jacobian matrix for h' at a point (0, m) has the form 

Dh'{0,m)= """^ ^ ° , 

^ ^ * A^^(777) 

with 

Nr{m) = Q"2(i(m - i)-i[r(47n - 3) + §(777 - i)] - Qim - i)-^) (30) 
= r{2m + 3){2m-l)-^Q-^ (31) 

where Q stands for Qr{l, fn) and where * denotes a term which depends on the 0(3) term 
in (j). In the case when m satisfies the fixed point equation ( l28ll this further simplifies 
(eliminating r) to give 

(2777 + 3)^ 

ISJm) = -. 

^ ^ 87772(2771-1) 

As the eigenvalues of Dh'{0, m) are A = (777 — |) and A' = Nj.{m) we easily conclude the 
following. 



Dynamics of an impact oscillator near a degenerate graze 18 

Proposition 4.4 The fixed point m^ with | < m^ < | has eigenvalues A, A' with 
< A < 1 and A' > 1. Also 

A-)-l, /i-)-l as r-)-l (ms -)■ |) 

while 

A— 7-0, yU— T-oo as r— )-0 (ms — )■ |). 
r/ie eigendirection for fi is the m-axis. D 

Note that the eigendirection for A would be the s-axis if there were no 0(3) terms in 0. 

Since (0,777.3) is hyperbohc it has a unique smooth stable manifold W^ transverse to 
the 777-axis. Then x{W^) is a smooth curve /3 through the origin in R^ and tangent to the 
s-axis from above (77 > 0) or below (77 < 0) according as a < or a > 0, and the same 
holds in the original (a, t) coordinates. By construction, the invariance of W^ under h' 
implies invariance of /3 under h, which in turn implies invariance of 7 := p{P)r\T\.Q under 
G. The dynamical properties of h' near the saddle point (0, 7773) carry over via x to the 
stated properties of G on and close to 7. This completes the proof of Theorem 14.11 D 

(h) The case a'ij^) < 

Here the horizon H for (r, t) 7^ (t*,0) = w^ locally lies in the region f < and so 
F : Hq — 7- V and hence also G : 11^ — )■ 11^ are continuous maps in a neighbourhood of 
w^ in Hq . In detail, we see that to first order F maps a neighbourhood T of w* in Hq 
onto a neighbourhood of the origin in the sector 

r : 3(T + t > , t > 

in R^ (coordinates on V ) which is then taken by (j) to the region 

T" : |a'(r,)(T2 <v<Q , a > 

in n^J": see Figure [6l Finally, R maps T" to a similar if narrower region RT" in 11^. 
From this geometry it is reasonable to suppose that under iteration of G all points of 
T' will tend to Ho . We now look more closely to verify that this is the case. 

Consider an arc a in V given in (a, t) coordinates by t = £0" with £ > and a > 0. 
From f lT5]l and f lTI?]) we have 

p(a, la) = {a, -k{U + f^)^^ + £aV(^, <y) (32) 

g(a, la) = (a(l + i), -rk{3£ + 2e^)a^ + ea^u{e, a)) (33) 

where /z, z/ are continuous (in fact smooth) functions for all i and sufficiently small a 
(depending on i). 

The images of a under p and q are locally the graphs of two smooth functions 
V = p{a) and v = q{a-). Choose 6 > and constants i^ > 0, L > such that 

\fi{ia)\<K \iy{i,a{l + e)-'^\<L 
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Figure 6. The actions oi p,(p and R close to (r, ,0) witli a'{T^.) < 0: licrc t = t^ + a. 
Tiie shaded area in the top left diagram represents the visible part of V, and the other 
shaded areas are its images under the relevant maps as indicated. 



for < i < —'^kr and < o" < 5. We then find that if 5 is small enough so that 
K5< -A;(3 + £)i(l-r) 



we have 



L5 < -rA;(3 + 2£)(l 



P{a) > -k{U + t)\{2 + r) 

q{a) < -rk{3e + 2f )(1 + £)"^|(1 + 2r) 



(34) 
(35) 

(36) 
(37) 



for < (7 < 6 (remember that A; < 0). Hence 

q{a)/p{a) < Mr(l + 2r)(2 + r)~^ 

where M = {3 + 2i){3 + £)~\l + £)-^. It is easy to check that M < 1 for £ > and also 
that r(l + 2r) < 2 + r for < r < 1. Hence for < a < 6 we have 

q{a) < Np{a) 

where < A^ < 1. Thus we obtain the expected result: 
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Proposition 4.5 Under iteration of G the arc p{a) tends C'-' -uniformly to the r-axis 
over the interval [0,5]. D 

Using the results from Section |3] we deduce : 

Corollary 4.6 For a G [0, 6] and t = t^ + a the stable manifold W^ extends backwards 
to a point (r', 0) G Ho with t' <t^ . D 

Hence we have the following local picture of the dynamics of the impact oscillator. 

Theorem 4.7 Let w^ = (r*, 0) G Ho 6e a tangency point with a'{T^) < 0. Then there is 
a neighbourhood U of w^, in Uq such that U^ := f/ fl 11^ is smoothly foliated by arcs of 
stable manifolds of points (r, 0) G Ho with r > r*, each arc having end-points on Uq on 
opposite sides of r* . See Figure \2\ 

Note in particular that each sufficiently small arc [wi,w*] C IIo is taken by G to a 
smooth arc e tangent to IIo at w^ in the region cr > 0, w > 0. 




Lq Wj 



W^ = (T^,0) 



Figure 7. Chatter close to w* = (r*, 0) with a'{Tf,) < 0. Points of Ho with t > Ti, are 
hniit points of nondegenerate chatter that becomes increasingly degenerate as r ^^ r* 
from above. The arc [wi,w»] on Ho is taken by G to the arc e tangent to Ho at w*. 



5. Further degeneracy: dynamics near a swan point 

Finally we consider the dynamics of G close to a swan point (r*, 0) ; here a{T^:) = a'{T^:) = 
but a"{n) ^ 0. 

As before, we work with coordinates r = t^ -\- a and we write {a, rj) = {s, sm) = 
x{s, m). Using flT7|) we see that the maps p and (p then become 

p : (s, m) ^ (s, £s^(6m + 4m^ + m^ + 0(s))) (38) 

: (s, m) h-> (s(l + m), -£s^(6m + Sm^ + Sm^ + 0{s))) (39) 

where C. = — i4a"(r*) and where the terms 0{s) denote smooth functions of (s,m) that 
vanish when s = 0. 
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In {s,m) coordinates the map h becomes h : {s,m) i— )■ (s, m) where 

s(l+m) =s (40) 

rs^{6m + 8m^ + 3m^ + 0{s)) = f{6m + Arfi^ + m^ + 0{s)) (41) 

and substituting the first equation into the second for s 7^ we obtain 

r(6m + 8m^ + 3m^) = (1 + mf{<orri + 4m^ + rrV') + 0{s). 

The graphs of the two real-valued functions 

/ : m h-^ r(6m + 8m^ + 3m^)(l + m)"^ (42) 

g : m ^^ 6m + Am? + m^ (43) 

intersect at a unique point m^, which lies in (—1,0). This proves the following. 

Proposition 5.1 The map h has a fixed point at (0,m*) and no other fixed point on 
the m-axis. 

It is easy to check from the graphs of /, g that f'{m^) > g'{m^), and also to verify that 

Dh{0,m^)= ^ % (44) 

* f{m*)/g'{m^) 

This implies the next result. 

Proposition 5.2 The map h has a hyperbolic saddle fixed point at (0,m*). The m-axis 
is the unstable manifold, while the stable manifold is a smooth curve W^ through (0,m^,) 
transverse to the m-axis. 

The interpretation of these results now depends on the sign of a"{T^). 

If a"(T^) < then a{T) has a local maximum at r* with a{T^) = 0, and so 
from [5]) we have ■^Vc{T,t) > for all (r, t) in some neighbourhood of w^ = (t*,0) 
in n except at w* . Therefore F is defined on a neighbourhood of w^, in 11 and is a 
smooth local homeomorphism onto a neighbourhood of w^ in V. Hence if /3 = x(W^*) 
then 7 = p{P) n Uq is a smooth invariant curve 7 for G terminating at the swan point. 
As is clear from (138|) this curve 7 has cubic tangency with IIq at w^. 

Theorem 5.3 Let w^, = (t^,,0) & Uq be a swan point with a"{T^) < 0. The map G 
has a unique invarant arc 7 with end-point at w^, , the arc being smooth and cubically 
tangent to the r-axis at w^, from the side t < t^^ and v > Q. The fixed point w^ is 
exponentially attracting for 6*17, and is weakly normally hyperbolic (repelling) for G. 
All points in a neighbourhood of w^ but not on 7 lie on stable manifolds of regular points 
o/IIo (nondegenerate chatter). See Figure\E D 

Thus with a"{T^) < all initial data close to w^, lead to nondegenerate chatter, with the 
exception of those on the curve 7 which exhibit a form of highly degenerate chatter. 

If on the other hand a"{T^:) > then we have Vc{T,t) < for sufficienatly small 
t > and T close to r*, and so the image under F of a neighbourhood of w* in U^ will 
include no points close to w^,. Therefore the local analysis of low- velocity impacts which 
is the subject of this paper will not apply. 
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Figure 8. Stable manifolds of chatter points close to a swan point (t*,0) where 
a'{u) = and a"{u) < 0. 



6. Unfolding the dynamics at the swan point 

Typically (the precise conditions on the functions / and 5^ in ([1]) are described in [16]) 
swan points will occur for only certain isolated values of c or of some other parameter 
in the system: swan points are codimension-1 phenomena. Therefore it is important to 
consider the bifurcations in the local dynamics as a generic transition through a swan 
point takes place. 

We shall take the parameter to be c for simplicity, although this is not crucial for 
the discussion. For clarity we also now re-introduce the c suffix in notation for the maps 
Pc, 4>c, Fc and Gc and the surface V/. We suppose that adr^) = a'^T^) = when c = Cq 
and let 






d{c, t) :- 

where a{c,T) = adr). As c increases through cq the function a^r) undergoes a Morse 
transition at which a pair of zeros ri, T2 of adr) is created or annihilated according to 
the sign of d^ := d{cQ,T^). Note that under the Correspondence Principle the creation 
of two tangency points at a swan point corresponds to the creation of two cusps at a 
swallowtail (Figure [2]^c)). 

Without loss of generality we suppose do < 0. Thus as c increases through cq an 
interval {ti,T2) is created on which ac{T) reverses its sign, and a new component of Pc 
(a swan configuration) is created out of the point w^ as described in [16] . See Figure [91 
The creation of the swan is in fact a version of the lips transition (Figure [21(a)) with 
particular geometry dictated by the context, a characteristic geometry of singular sets 
noted earlier by Bruce [llj . 

There are two cases to consider, according to whether o-o (''"*) '■— '^co(t*) > {positive 
swan point) or ao(r*) < {negative swan point). 
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Figure 9. Graphs of a{T) and associated swan configurations in the plane 11 in two 
cases: (a) positive, (b) negative. 





T<X, 



T<T 



T=T, 



Figure 10. Sections through V^ in the {t, w)-plane for increasing values of r with 

T < T^. 



6.1. Positive swan point: a"{T,,) > 

For c < Co there are no low- velocity impacts occuring close to w*, and for c = Cq although 
the acceleration vanishes at r^, there are still no low- velocity impacts. We therefore focus 
on the low- velocity dynamical phenomena that are created when c > Cq . Once again 
we now suppose such c fixed and dispense with it as a suffix. 

To help picture the dynamics of G (and in particular the geometry of F) we show in 
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Figure 11. Positive swan configuration: the visible part 5* of V is bounded by an arc 
of the horizon H and an arc of the zero locus Z. See text for details. 



Figure [TO] a sequence of sections through V at constant r for various choices of r < r*. 
This information can be obtained from the Taylor expansion (ITTl) . regarding Vc{T^, + a, t) 
as a function of t unfolded by the parameter a. Analogous sketches for r > t^, are 
obtained by reversing the signs of both t and v. 

In Figure [11] we show the horizon H and the zero set Z in (r, t)-coordinates on 
V , thus viewed from the 'vertical' v direction. Here A = (ti,0) and B = (t2,0) where 
Ti < T2 are the two zeros of a^r) near r* for c > cq and close to cq. The shaded area S 
indicates the connected component of the visible part of V that contains the interval 
AB (see Section [22]) . 



As already seen in Chapter [H the dynamics of G close to A and B are as in Figure [5] 
and Figure [7| respectively. These combine to give the swan dynamics in the way that 
we now describe. 

Figure [T2t^ a) shows the same configuration as Figure [TT] but viewed along the t- 
axis, or in other words its image in 11 under the projection p . Points p{C),p{D), . . . are 
denoted C, D', . . . ; note that A' = A and B' = B. The arc BE of H is taken by p to 
a smooth arc BE' in 11^ tangent to Ho at B and transverse to Ho at E': together with 
the interval E'B of Ho it bounds the simply-connected region S' = p{S) C Hq. 

Figure [T2l^ b) shows the image of the configuration in Figure [11] under the re-set 
map 0, with points (/>(C), 0(-D), . . . denoted C", D", .... The arc AE of Z is taken by (j) 
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to a smooth arc AE" in n^j" tangent to Ho at A and transverse to Ho at E" (part of the 
apparent outhne in view of the Correspondence Principle: see Proposition 12. 3p : together 
with the interval AE" of Ho it bounds the simply-connected region S" = 0(5) C 11^ . 

To construct the dynamics of G we compose the inverse of p with and then with 
the restitution R. The latter takes the arc AE" in Hq to a smooth arc in Uq tangent 
to Ho at A and transverse to Ho at E". Its intersection with S' consists of a smooth arc 
e also tangent to Ho at A and with its other end-point N lying on the arc LB of the 
apparent outline. 

The G-invariant curve 7 at the degenerate chatter point B is tangent to Ho at B 
and transverse to Ho at its other point of intersection. It separates S' into two regions: 
one foliated by stable manifolds of points of AB (nondegenerate chatter), the other 
consisting of points that eventually leave S'. See Figure | 



6.2. Negative swan point: ao(r*) < 

For c < Co we have a{T) < for all r in a neighbourhood of r* and so Theorem 13.11 
applies: every point lies in the stable manifold of some point of Ho and nondegenerate 
chatter takes place locally as in Figure |3]^b). 

For c = Co Theorem 15.31 applies: there is a smooth stable manifold for (r=K,0) that 
is cubically tangent to Hq and corresponds to highly degenerate chatter, while all other 
orbits experience nondegenerate chatter becoming increasingly degenerate as r — )■ r^, 
from either side as in Figure [HI We therefore focus on the case c > Co . 

The analogue of Figure [10] for this case is obtained by inverting the t^-axis. In 
Figure [Hj we again show the horizon H and the zero set Z as viewed from the v- 
direction, with A,B as before corresponding to the values ti,T2 of r where ai^r) = 0. 
Now, however, it is a different region T = 7i U 72 U Ts that is visible, bounded by the 
arc BC of Z , the arc DA of H , and the (dotted) arc DC which is the shadow of DA, 
that is the set of "next-hit" points starting from DA: more precisely it is the set of 
points (r, v; ti) G V for which there exists tg > with (r, v; to) £ DA but (r, v; t) ^ V 
for all tQ < t < ti. The visible part of V is the region lying outside the loop ADCBA 
(including AD but not including the open arc DC). 

Under the correspondence principle in Figure [HJ the points A, B each correspond 
to cusp points elsewhere on H: in fact B corresponds to D while A corresponds to 
its counterpart (not shown) with t < 0. Figure [T5](a) shows the images A', . . . , £" of 
the points A, . . . ,E of V under projection p, while Figure [T^ b) shows their images 
A", . . . , E" under the re-set map . The only visible part of H is the arc from A to 
D: under p this projects to the arc AD' in IIq while under it is taken to the interval 
A"D" on the r-axis Ho. 

At points of the shadow arc CD apart from D the velocity x^ is nonzero (in fact 
it has to be negative) and so takes CD to an arc C"D" with f < except at D"; 
moreover C" is a cusp point of P by Proposition 12.21 From this we see the nature of 
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Figure 12. The images under (a) the projection p and (b) the re-set map of various 
regions of V' in Figure [TT] for an unfolded positive swan point: see text for details. 



the discontinutity of the map (j) o F (and hence of G): cutting n"*" along the arc AD' 
then mapping one edge (heavy) to the interval A"D" of Ho and the other edge (dotted) 
to the arc D"C" . Meanwhile the interval AB of Figure [T5l(a) is mapped by F to the arc 
CB of Z which is then taken by to the arc C"B of Figure ITBT b). 

These pictures repay contemplation in order to form a clear mental picture of the 
dynamics of G . Figure [16] shows schematically how (poF takes each segment of the swan 
region into the next segment; remember, however, that the map G also incorporates the 
restitution R at each step and takes each of 11^, Hq into itself. We now describe the 
dynamics. 

In Section 14.11 we have already obtained local models for the dynamics of G close 
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Figure 13. Invariant foliation for dynamics of G in the region S' . Outside S' the 
unfolding of the swan point gives no local information. 
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Figure 14. Negative swan configuration on V: the visible part (shaded) of V here 
is bounded by an arc of H and an arc of Z together with an arc DC of the shadow of 
H. 



to the two tangency points A and B (Figure [5] and Figure [7] respectively). As there are 
no discontinuities of G close to B, the neighbourhood of B that is foliated by stable 
manifolds W^ of points (r, 0) with t > T2 (by Theorem 14. 7p must extend at least as far 
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(b) 




Figure 15. The images under (a) the projection p and (b) the re-set map of various 
regions of V' for an unfolded negative swan point: see text for details. 



as the stable manifold W^,^ = W^*(X) that contains the point A, where X = {t^, 0) with 
^3 > T2. See Figure ITT! 

Proposition 6.1 The stable manifold W^{X) is tangent to Uq at A. 

Proof. We argue that F{W^{X)) C V is tangent to Z at C, from which the result 
follows by applying p . 

Since the map has a cusp singularity at C by Proposition I2.2[ if F{W^{X)) were 
not tangent to Z at C then o F{W{X)) would be aligned with the direction of the 
cusp at C" in Figure [H In that case i? o o F{W'{X)) = G{W'{X)) = W'{X) would 
be tangent to the arc e' = R{D"C") at its endpoint R{C"), and this is not consistent 
with the geometry of Figure [3 stable manifolds of points of Ho to the right of and close 
to B are transverse to e'. 
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Figure 16. Schematic representation of tlie action oi (f) o F within the swan region. 
Matching arrows indicate arcs and their images. The third step is discontinuous along 
AB. 



D 

Let A denote the arc of W^{X) from A to Q = W''^{X) fl e'. The inverse image 
G~^{\) is a closed loop around the discontinuity arc 6 = AD' and with both ends at A. 
The inverse images of S under iterates of G form a sequence of arcs tangent to Ho at 
the degenerate chatter point A and accumulating on its invariant arc 7, and therefore 
the inverse images of A form a sequence of loops also accumulating on 7. 

Stable manifolds W^(Y) of points Y to the right of X on Ho (i.e. with r > rs) 
experience no discontinuities locally: those for Y close to X oscillate more and more 
wildly under backwards iteration of G as they attempt to adhere to the loops of G'"'{X). 
See Figure [T71 



7. Global dynamics 

The dynamical phenomena that become apparent in the unfolding of a negative swan 
point give some clues about the global geometry for dynamics of an impact oscillator. 
In particular we see that the particular geometry of the smooth map p o 0~^ close 
to a tangency point w* = (r*,0) G Ho (so that z = 0~^(w*) is a cusp point of H) 
plays a key role. At a tangency point w* where a'{T^) > the nearby stable manifolds 
W^ with T > T^: have infinitely many loops accumulating on the invariant curve 7, 
that is the stable manifold of w* Therefore we see that the global behaviour of the 
stable manifolds of degenerate chatter points, and in particular their intersections with 
unstable manifolds of fixed or periodic points elsewhere in the phase space n^j" , are 
crucial in organising the global effects of the complicated dynamics unleashed at a 
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Figure 17. Dynamics of G created at unfolding of a negative swan point: see text for 
details. 



(negative) swan bifurcation. The swan configuration created at a negative swan point 
will persist until variation in the parameter c causes it to impinge on some other branch 
of the apparent outline P . Careful study of the generic ways in which this can happen 
will provide important tools for understanding the complicated global structure of 
discontinuous planar discrete dynamical systems modelling 1-degree of freedom impact 
oscillators. 
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